





PURIFICATION OF IMPURE DENSITY OPERATORS AND THE RECOVERY 

OF ENTANGLEMENTS 

V. I. Man'koQ, G. Marmo, E. C. G. SudarshanQ and F. Zaccaria 

Dipartimento di Scienze Fisiche 



O I Universita di Napoli "Federico II" and 

00 . 

Istituto Nazionale di Fisica Nucleare, Sezione di Napoli 

> '. Mostra d'Oltremare, Pad. 20 - 80125 Napoli, Italy 

00 
O, 

O , Abstract 



is determined in terms of phase dependent multiplication. 



The need to retain the relative phases in quantum mechanics implies an addition law 
parametrized by a phase of two density operators required for the purification of a density 
§ : "^*""- ^'" " ^'°"" "^'' ^"^"*"" tomography and the Wigner function. Entanglement 

> 

^: 

c^ ■ Quantum mechanics as traditionally formulated [1] involves three principles. The states of the 

system are normalized vectors in a Hilbert space, (selfadjoint) linear operators correspond to (real) 
dynamical variables; and the expectation value for any dynamical variable is bilinear (sesquilinear) 
in the state vectors. But the overall phase of the state vector is irrelevant in those computations 
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involving only that one state. The state corresponds to a "ray" in Hilbert space {ip : e^^ipo} ■ 
Alternatively the expectation value can be expressed by means of the density matrix 

which is defined on the ray. 

These density operators are "pure" and of rank one: 

p^ = p; p>0; trp = l; p"^ = p ■ 

There is a probability addition of the density operators. If pi and p2 are two pure density matrices 

p = cos^ 0pi + sin^ Qp2 

also is a density matrix; but it is not pure, but mixed. These are the interior points of the convex 
set of density operators. Hence 
So 

P^T^P, p-p^>0. 

Can we have an additional operation of constructing a pure density operator from the sum of the 
two density operators pi, p2? We must have such a construction since the superposition principle 
of quantum mechanics [1] tells us how to add two state vectors ipi and 1^2 to form a superposed 
state vector 

V' = I cos Q ^pi + e**^ sin Q ^2] 

which is normalized if ipi, 1^2 are orthonormal. For nonorthogonal states, we have to use the 
normalization factor 

(1 + sin 29 cosv? | (^1 1^2) |)"^^^ 



where y? includes the phase of {ipi \ 1^2) along with ip^. In terms of density operators, the 99-addition 
law which we introduce is 

p{x, y,(p) = {cos Q^pi (x) + e*''' sin 6-?/'2 (x) } {cos Q^pl (y) + e**^ sin Q^pl in) } 

= cos^ Qpi{x, y) + sin^ Qp2{x, y) + sin26 cos(ppi2{x, y, ip). 

In terms of the Wigner function[2], we get a (/9-addition law which has an interference term 
proportional to sin 20 cosy?-, i.e. 

^ {I1 p) = cos^ Wi (g, p) + sin^ W2 {q, p) + sin 20 cos ip I12 {q, p, ^p) 

where I12 is the (generalized) Wigner function corresponding to the operator having the struc- 
ture of the root square of a convolution of the product of the two density matrices, which we 
symbolically denote as 



There is thus a one-parameter addition law of density operators and of Wigner functions with 
probabilities cos^ and sin^ and with extra intereference term. Note that this W (g, p) is pure 
and satisfies the purity criterion 



j jdpdq{W{q,p)Y = ^. 



For quantum tomography [3], also we can construct superposition of tomograms using a one 
parameter addition law: 

$ (A, /x; x) = cos^ $1 (A, p; x) + sin^ $2 (A, p; x) + sin 20 cos ip $12 (A, p; x) . 

Here tomograms $i.$2 determine the probability density of quadrature x = Xq + pp in the 
pure states ipi and 1^2- The tomogram $12 corresponds to the interference term J12. 



The passage from the impure density operator cos^ Opi + sin^ Op2 to the pure (^-dependent 
addition may be called purification. Note that the purification introduces the relative phase (p 
which was not in pi or p2- 

The density operator of a composite system AB with subsystems A and B may be chosen pure 
or impure. For a pure density operator p^s, one can get the density operator p^ and ps by the 
partial trace operation 

Pa = trij(p^ij); PB = tr^(p^B). 

It is not necessary that A and B have the same dimensionality. Unless pab is a direct product of 
pure states of A and 5, a pure pab yelds impure pA and ps- But they will have the same rank R 
and the same nonnegative eigenvalues which sum up to unity. The density operator 

p'ab = Pa® Pb^ Pab 

is impure. Thus, the whole is greater than the parts: there is additional information in Pab- These 
are the "entanglement" terms[4]. 

We purify the product p^^ in the same way that we used before for the mixture of two density 
operators pi and p2. Here we have n such pure states mixed together and need (^ — 1) phase angles 
(fi = 0, (f2, ■ ■ ■ V^ra- The diagonal form of p'ab has only n nonzero diagonal elements. We need to 
introduce the offdiagonal elements 

in the (j, k) location. Note that while we have n{n — l)/2 offdiagonal terms, there are only (n — 1) 
phases (fj. 

The purification of the density matrix p^^ we call as the ip- multiplication law of the density 
matrices pA and p^. 



While purification of an impure density addition is dependent on one pliase angle, the form 
of the entanglement is constructed depending on i? — 1 phase angles. These have to be obtained 
from other considerations. 

The same kind of (^-addition law and (^^-multiplication law holds for other representatives like 
the quantum tomograms, the diagonal coherent state distribution function in quantum optics [5] 
and the Husimi - Kano[6] density of coherent state projection operators. We expect to return to 
this discussion elsewhere. 
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